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ABSTRACT 


While  the  mechanisms  of  flow  separation  for  a  submarine  in  steady  rectilinear 
flow  are  well  understood,  the  flow  separation  behaviour  of  a  submarine  during 
a  manoeuvre  is  not  as  well  understood.  The  importance  of  this  has  led  us  to 
investigating  the  flow  separation  over  a  submarine- like  body  of  revolution  in  a 
turn  of  fixed  radius. 

The  flow  separation  over  the  manoeuvring  body  is  examined  by  using  a 
deformed  body  placed  in  rectilinear  flow  to  simulate  the  flowfield  created  by  a 
body  moving  through  a  circular  path.  This  method  was  used  by  Von  Karrnan 
at  the  Guggenheim  Airship  Laboratory  and  by  Gurzhienko  in  Moscow  in  the 
1930s  to  investigate  the  manoeuvring  of  airships. 

The  deformed  model  is  constructed  such  that  the  local  angle  of  attack  is 
the  same  for  a  regular  model  placed  in  curved  streamlines.  Velocity  should 
be  varied  linearly  across  the  test  section  to  ensure  the  local  velocities  at  the 
deformed  model  surface  are  equal  to  the  velocity  along  the  manoeuvring  model. 

Numerical  simulations  conducted  using  the  CFD  code  FLUENT  show  the 
validity  of  this  transformation.  Results  for  a  body  placed  in  a  rotating  reference 
frame  are  compared  to  a  deformed  body  in  rectilinear  flow. 
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Flow  Over  a  Body  of  Revolution  in  a  Steady  Turn 


EXECUTIVE  SUMMARY 

The  problem  of  obtaining  detailed  data  of  the  flow  over  a  body  of  revolution  in  a 
steady  turn  is  approached  in  an  innovative  way.  A  co-ordinate  transformation  is  applied 
to  produce  a  curved  model  which,  when  placed  in  a  steady  rectilinear  flow,  re-creates  the 
flow  geometry  of  a  straight  body  in  curved  flow.  Since  the  curved  model  is  sitting  in  a 
steady  flow  situation  (which  can  be  created  with  a  wind  or  water  tunnel)  detailed  flow 
measurements  of  the  vortices  created  by  turning  objects  can  be  performed. 

In  this  report,  an  analysis  of  the  co-ordinate  transformation  is  presented.  Preliminary 
numerical  studies  have  been  conducted  with  the  CFD  code  FLUENT.  A  body  of  revolution 
was  placed  in  a  rotating  reference  frame  to  simulate  the  flow  over  this  body  in  a  gentle 
turn.  The  results  of  this  simulation  were  compared  to  the  corresponding  curved  body 
(whose  geometry  is  related  to  the  turn  radius  and  angle  of  attack  of  the  turning  body)  in 
rectilinear  flow.  Results  were  obtained  for  both  the  surface  flow  features  (such  as  average 
static  pressure  and  skin  friction  distributions)  and  downstream  flow  features  (the  formation 
of  the  main  off-body  vortex). 

The  results  show  that  the  distributions  of  skin  friction  along  the  curved  and  straight 
models  differ  slightly,  in  that  we  observe  more  flow  separation  for  our  curved  bodies. 
This  gives  a  correspondingly  lower  average  pressure  for  the  curved  body.  In  addition,  the 
strength  of  the  vortex  produced  by  the  curved  body  is  lower.  Further  numerical  work  will 
investigate  the  effects  of  the  coriolis  and  centrifugal  forces  on  the  formation  of  the  off-body 
vortex. 

These  results  show  at  this  stage  that  any  experiments  incorporating  this  method  of 
curved  bodies  will  be  able  to  replicate  the  dominant  flow  features  that  occur  over  a  straight 
body  performing  a  steady  turn. 
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1  Overview  and  Literature  Survey 

1.1  Introduction 

The  behaviour  of  a  fluid  flowing  over  a  body  in  a  manoeuvre  is  a  problem  that  is 
present  in  many  fields  of  research. 

Most  of  the  research  into  the  flow  fields  created  by  a  manoeuvring  object  has  been 
related  to  the  calculation  of  dynamic  stability  derivatives.  These  derivatives  are  required 
during  the  design  process  and  performance  evaluation  of  any  ob  ject  designed  to  execute 
a  controlled  manoeuvre  in  a  fluid,  be  they  aeronautical  in  nature  (aircraft,  missiles,  bal¬ 
listic  shells)  or  maritime  (ships,  submarines,  torpedoes).  When  the  object  in  question  is 
considered  to  execute  a  turn  with  a  fixed  radius,  the  results  can  be  applied  to  fields  such 
as  turbomachinery  and  wind-turbines. 

A  greater  understanding  of  unsteady  fluid  flow  in  recent  years  has  focused  research 
into  the  forces  on  the  structure  of  the  body  caused  by  unsteady  fluid  flow.  In  particular, 
studies  of  vortex  breakdown  have  been  carried  out  to  investigate  the  induced  vibration  of 
aircraft  control  surfaces. 

Unsteady  flow  becomes  a  major  problem  for  maritime  craft  as  it  can  induce  vibrations 
of  the  propeller,  caused  by  the  sudden  changes  in  the  effective  angle  of  incidence  of  the 
propeller  blades  as  they  rotate  through  the  unsteady  flow.  Vortices  themselves  create 
excessive  hydro-acoustic  noise.  In  the  specific  case  of  a  submarine’s  acoustic  signature,  it 
is  of  great  importance  to  minimise  and  preferably  eliminate  any  flow  separation  in  order  to 
reduce  the  hydro-acoustic  noise  from  the  overall  acoustic  signature.  While  the  mechanisms 
of  flow  separation  are  well  understood  for  a  submarine  in  steady  flow,  the  flow  separation 
behaviour  (and  hence  the  changes  in  the  hydro-acoustic  signature)  of  a  submarine  during 
a  manoeuvre  are  not  as  well  understood. 


A 


D 


Figure  1:  System  of  vortices  produced  by  a  submarine  in  a  manoeuvre  (taken  from  Lloyd 
and  Campbell) 


Figure  1  shows  the  system  of  vortices  produced  by  a  submarine  during  a  turn.  Each 
appendage  on  a  submarine  will  produce  its  own  system  of  vortices,  such  as  the  sail  (B), 
tailplane  (C)  and  rudder  and  tail-fins  (D).  The  strength  of  these  vortices  is  far  less  than 
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those  cast  off  the  body  itself  (A).  These  two  vortices  are  produced  by  the  cross-flow 
separation  over  the  body  of  revolution.  As  the  angle  of  attack  increases  along  the  body 
(see  figure  4),  the  fluid  flows  over  the  outboard  side  of  the  body  to  the  inboard  side.  At 
this  point,  the  flow  is  no  longer  aligned  with  the  body  of  revolution,  but  is  in  fact  flowing 
across  the  circular  hull  cross-sections.  At  some  point  along  the  body,  this  flow  will  start 
to  separate,  and  an  off-body  vortex  is  produced. 

1.2  Analysis  of  a  rudder-induced  manoeuvre 

Lewis  [1]  provided  an  extensive  review  of  controllability  and  model  testing  in  relation 
to  issues  of  Naval  Architecture.  Chapter  Nine  of  this  volume,  compiled  by  C.L.  Crane,  H. 
Eda  and  A.  Landsburg  contains  a  detailed  study  of  the  analysis  of  a  ship  in  a  turn. 


Figure  2:  Path  of  an  object  in  a  rudder-induced  turn  (taken  from  Lewis) 

Consider  a  ship  moving  through  the  water  in  a  straight  path  with  a  steady  velocity 
(figure  2).  The  ship  experiences  forces  only  along  its  :r-axis.  The  thrust  force  produced  by 
the  propeller  is  balanced  by  the  viscous  drag  forces  acting  on  the  hull  and  its  appendages. 
When  the  rudder  is  deflected,  it  produces  a  force  which  has  a  component  along  the  y- axis 
of  the  ship.  For  a  starboard  turn,  the  rudder  is  deflected  towards  starboard,  producing  a 
force  in  the  port  direction,  pointing  away  from  the  centre  of  turn.  This  force  produces  a 
moment  that  causes  the  ship  to  rotate  about  its  2-axis.  The  stern  of  the  ship  will  move  to 
port,  while  the  bow  of  the  ship  will  rotate  towards  the  centre  of  the  turn.  This  rotation 
causes  the  bow  of  the  ship  to  point  into  the  turn  with  a  positive  angle  of  attack,  or  drift 
angle  /J.  The  hydrodynamic  forces  acting  on  the  ship’s  hull  which  is  placed  at  a  positive 
angle  /3  to  the  flow  will  produce  a  force  which  is  directed  towards  the  centre  of  turn.  This 
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force  will  create  a  change  in  the  lateral  movement  of  the  center  of  gravity.  In  general,  as 
the  drift  angle  ft  increases,  the  hydrodynamic  force  acting  on  the  hull  towards  the  turn 
centre  will  increase  as  well.  Therefore,  in  order  to  perform  a  tighter  turn,  the  drift  angle 
must  be  large. 

When  the  rudder  remains  in  a  fixed  position,  a  steady  turn  will  develop.  In  this 
situation,  the  hydrodynamic  and  inertial  forces  and  moments  are  balanced,  and  the  path 
of  the  centre  of  gravity  is  a  circular  arc. 

The  immediate  response  of  a  ship  to  the  deflection  of  the  rudder  produces  a  transient 
turning  profile  where  large  surge,  sway  and  yaw  accelerations  occur.  Since  this  project  is 
concerned  with  the  fluid  flow  over  a  body  in  a  steady  turn,  the  initial  transient  behaviour 
will  not  be  considered. 

The  location  of  the  pivot  point  in  a  steady  turn  locates  the  position  on  a  ship  where 
the  flow  of  water  past  the  ship  is  parallel  to  the  ship’s  x-axis.  The  location  between  the 
pivot  point  and  the  centre  of  gravity  can  be  given  by  xc  =  Rsinfi,  where  R  is  the  turn 
radius.  Since  small  radius  turns  are  usually  associated  with  large  drift  angles,  and  large 
radius  turns  with  small  drift  angles,  the  product  Rsinft  does  not  vary  significantly  for 
different  ships,  or  indeed  the  same  ship  performing  turns  of  different  radius.  For  most 
ships  the  pivot  point  lies  between  the  bow  and  1/5  aft  of  the  bow  (see  Mandel  [2]).  Based 
on  empirical  data,  the  drift  angle  ft  can  be  computed  either  as: 

ft  =  22.5L/R  +  1.45  (1) 

or 

ft  =  18  L/R  (2) 

Using  equation  1  gives  values  of  xc  ranging  from  0.4  to  0.5L.  Equation  2  gives  xc  =  0.3 L. 

Once  the  ship  has  settled  into  the  steady  part  of  the  turn,  the  accelerations  in  the 
radial  direction  r  and  along  the  y-axis  v  are  zero.  The  linearised  equations  of  motions  for 
a  steady  turn  can  be  written  as  : 


-YVV  -  (Yr  -  An,)  =  Y65r 

(3) 

-Nvv  -  Nrr  =  Ns6r 

(4) 

where  : 


Y  =  force  along  the  y  axis  of  the  ship  (sway) 

N  =  moment  about  the  z  axis  of  the  ship  (yaw) 
v  =  velocity  along  the  y  axis  of  the  ship 
r  =  rate  of  rotation  about  the  z  axis  of  the  ship 
u\  =  initial  equilibrium  value  of  x-component  velocity 


Yr 


-  §y 

dr 


Yv 


—  §y 

dv 
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N  —  — 

iV v  “  dv 
A  =  mass 


8r  =  rudder  deflection  angle 


If  the  stability  derivatives  Yr ,  Yv ,  ATr,  A^  and  control  derivatives  and  A^  are  known, 
then  equations  3  and  4  can  be  solved  for  r  and  v.  The  solutions  show  that  the  steady 
turning  radius  is  proportional  to  ship  length  and  inversely  proportional  to  8r.  and  that  /3 
is  directly  proportional  to  8p.  This  linear  analysis  is  useful  for  estimating  steady  turning 
radii  and  /?  for  large  diameter  turns,  but  for  ships  in  a  tighter  turn,  the  validity  of  the 
linear  equations  of  motion  breaks  down. 


i  z  axis 

1  z  axis 

C.G. 

Av+Ax  r 
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-=■ 

—  y  axis 

-=■ 

-= 

Ys5 

Y.v  +  Y  r 

V  r 

Y  Vv  +  Yr  r 
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0  K 

J 

) 

■=-  y  axis 


(A)  At  initial  rudder  deflection  (B)  At  steady  turning  radius 

Figure  3:  Location  of  forces  for  a  ship  in  a  starboard  turn  looking  aft  (taken  from  Lewis) 


The  use  of  a  rudder  also  introduces  pitching  and  rolling  motion  (see  figure  3).  Note 
that  part(B)  of  figure  3  is  a  graphical  representation  of  equation  3.  The  rolling  motion  in  a 
steady  turn  is  created  by  the  moment  caused  by  the  rudder  force  {YsSr)  which  points  away 
from  the  turn  centre,  and  the  force  on  the  hull  ( Yvv  +  Yrr)  acting  towards  the  direction 
of  turn.  A  ship  will  tend  to  have  outboard  roll  angle  during  a  steady  turn  (i.e.  lean  away 
from  the  turn),  whereas  a  submerged  submarine  will  always  lean  into  the  turn.  This  is 
due  to  the  location  of  the  line  of  the  hull  force  Yvv  +  Yrr,  which  because  of  the  presence  of 
the  sail  is  considerably  higher  relative  to  the  centre  of  gravity  than  for  a  surface  ship.  The 
sail  on  a  submarine  acts  as  a  lifting  device  during  a  steady  turn,  so  the  sail  contributes 
to  the  magnitude  and  height  of  the  Yvv  force  which  increases  the  rolling  moment.  Due  to 
this  increase  of  Yvv  force  a  submerged  turning  submarine  will  initially  experience  a  large 
inboard  roll  called  a  snap  roZZ,  followed  by  smaller  inboard  rolls.  Arentzen  and  Mandel 
[3]  investigated  the  importance  of  the  sail  configuration  in  the  amount  of  roll  induced  in 
a  turning  manoeuvre,  concluding  that  the  presence  of  a  sail  increases  the  roll  excitation 
in  a  turn,  but  also  increases  the  magnitude  of  roll  damping,  and  therefore  dampens  the 
amplitude  of  the  overshoot  of  the  snap  roll. 
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1.3  Manoeuvring  characteristics  of  bodies  of  revolution 


There  exists  a  large  body  of  work  in  the  open  literature  regarding  the  nonlinear  aero¬ 
dynamic  performance  of  bodies  of  revolution.  In  particular  this  work  has  been  approached 
from  the  area  of  ballistic  research,  and  is  mainly  concentrated  on  the  calculation  of  pitching 
and  rolling-moment  calculations  for  ballistic  projectiles.  A  lot  of  theoretical  and  experi¬ 
mental  work  was  conducted  at  NASA  Ames  Research  Centre  in  the  late  1960s  and  early 
1970s  (see  Tobak  et.  al  [4],  Seliiff  and  Tobak  [5],  Levy  and  Tobak  [6]).  This  experimen¬ 
tal  work  focused  on  the  measurement  of  forces  and  moments  for  bodies  of  revolution  at 
high  Mach  numbers.  Flow-visualisation  work  was  also  performed  using  the  vapour-screen 
method,  revealing  the  dependence  of  some  moment-coefficients  on  the  appearance  of  dis¬ 
placed  vortices.  Schiff  [7]  used  a  simple  numerical  scheme  to  solve  for  inviscid  flow  of  a 
perfect  gas  around  a  body  performing  both  pitch  and  roll  in  supersonic  flow.  Schiff  used 
a  rotating  co-ordinate  frame  to  simulate  the  manoeuvring  of  the  projectile.  The  flow  was 
steady  within  the  rotating  frame,  allowing  the  steady  Euler  equations  to  be  solved.  Good 
agreement  was  achieved  with  experimental  data  for  moderate  values  of  crossflow  flight  ve¬ 
locity.  Schiff  s  work  was  carried  further  by  Weinacht  and  Sturek  [8]  and  Weinacht  [9]  who 
used  the  thin-layer  Navier-Stokes  equations  to  compute  pitch  damping  coefficients  for  par¬ 
ticular  projectile  geometries.  Further  work  using  rotating  reference  frames  was  conducted 
by  Buckland  [10]  and  Lesage  [11].  Lesage  used  the  commercial  CFD  code,  TASCflow,  to 
solve  the  Reynolds-averaged  Navier-Stokes  equations,  using  the  k  —  e  turbulence  model 
and  wall  functions  to  calculate  the  Reynolds  stresses  and  thermal  diffusion  for  turbulent 
solutions. 


1.4  Manoeuvring  characteristics  of  ships  and  submersibles 


Extensive  research  has  been  conducted  in  investigating  the  manoeuvring  characteristics 
of  ships,  and  more  recently,  submersibles.  A  lot  of  this  research  has  focused  on  the 
methods  of  ship  design,  and  in  particular  methods  to  predict  ship  performance  at  an 
early  design  stage.  Work  conducted  by  Burcher  [12,13]  describes  the  various  methods  of 
model  testing  used  to  measure  stability  and  control  characteristics  of  ship  forms.  Burcher 
describes  free  model  tests  using  radio  controlled  scale  models,  and  constrained  model  tests. 
Constrained  model  tests  included  both  rotating  arm  tests  and  oscillatory  testing  using  a 
Planar  Motion  Mechanism.  Fujino  [14]  examined  the  effects  of  frequency  dependence  on 
stability  derivatives  for  ship  motion. 

Feldman  [15]  conducted  rotating  arm  and  oscillatory  experiments  to  investigate  stabil¬ 
ity  and  control  characteristics  of  submarines  and  submersible  vehicles.  Mendenhall  et.  al 
[16]  conducted  a  thorough  investigation  into  the  vortex  shedding  from  circular  and  non¬ 
circular  bodies  at  high  angles  of  attack.  The  analytical  methods  he  developed  to  predict 
vortex  shedding  were  used  to  predict  both  the  aerodynamic  characteristics  of  missiles,  and 
the  hydrodynamic  characteristics  of  submersibles  in  a  manoeuvre  (see  Mendenhall  and 
Perkins  [17]). 
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1.5  Vorticity  measurements  of  a  manoeuvring  submarine 

A  thorough  experimental  program  has  been  conducted  by  Lloyd  [18]  and  Lloyd  and 
Campbell  [19].  Lloyd  proposed  the  use  of  a  simple  vortex-based  computer  program  SUB- 
SIM  which  would  provide  an  early  approximation  of  the  forces  and  moments  on  a  sub¬ 
marine  hull  during  a  manoeuvre.  The  aim  of  SUBSIM  was  to  bridge  the  gap  between 
the  traditional  stability  derivative  approach  to  the  prediction  of  manoeuvring  characteris¬ 
tics,  and  the  use  of  CFD  (which  at  that  time  was  not  mature)  to  estimate  the  forces  and 
moments  on  a  body. 

The  flow  over  a  manoeuvring  submarine  is  dominated  by  vortices  which  are  shed  from 
the  hull  and  appendages  (such  as  bow  planes  and  the  sail).  The  strength  and  position 
of  the  appendage  vortices  can  be  predicted  using  lifting  line  theory.  The  body  vortices 
are  affected  by  the  angle  of  incidence  and  rate  of  turn  of  the  submarine.  SUBSIM  was 
designed  to  use  empirical  formulae  to  represent  the  positions  and  strengths  of  these  body 
vortices. 

Lloyd  tested  three  submarine-like  1  bodies  of  revolution  in  a  water  tunnel  in  both 
rectilinear  and  rotating  flow.  He  used  a  probe  developed  by  Freestone  [20]  to  measure 
the  strength  of  the  vortices  cast  from  the  hull.  Measurements  of  vorticity  were  taken  at 
10  degree  angular  intervals  and  25  mm  radial  steps  from  the  body  surface.  Experimental 
variables  such  as  angle  of  incidence  and  turn  radius  were  varied.  By  varying  the  angle  of 
incidence,  Lloyd  was  able  to  shift  the  pivot  point  along  the  model,  although  for  the  models 
in  a  tight  turn,  he  still  had  large  local  angles  of  incidence  at  the  nose.  This  contradicts  the 
work  of  Mandel  [2]  who  showed  that  the  local  angle  of  incidence  at  the  nose  approached 
zero  as  the  turn  radius  decreased.  Lloyd’s  experiments  showed  that  for  a  body  executing 
a  turn  with  a  large  radius,  the  flow  approached  from  the  outboard  side  of  the  body,  and 
vortices  were  shed  from  the  inboard  side.  When  the  body  was  placed  in  a  tighter  turn,  and 
the  pivot  point  located  aft  of  the  bow,  the  local  flow  direction  at  the  nose  was  from  the 
inboard  side  of  the  body,  and  over  the  tail  sections  of  the  hull,  the  flow  direction  was  from 
the  outboard  side.  Lloyd  hypothesised  that  this  would  produce  areas  of  mixed  vorticity, 
with  the  vortices  cast  from  the  nose  convecting  along  the  body,  and  then  interacting  with 
the  vortices  of  opposite  sign  which  were  cast  from  the  tail.  His  results  show  some  evidence 
of  this.  Lloyd  concluded  that  the  flow  curvature  generally  increases  the  vortex  strengths 
but  did  not  noticeably  shift  the  location  of  the  vortices.  The  presence  of  a  secondary 
system  of  vortices  of  opposite  sense  generally  located  on  the  outboard  side  of  the  body 
were  responsible  for  limiting  the  total  amount  of  circulation  in  tight  turns. 

Lloyd  held  the  body  in  a  fixed  attitude  while  changing  the  turn  radius.  It  should  be 
noted  that  in  a  tighter  turn  of  smaller  radius,  the  centrifugal  forces  will  be  larger  (assuming 
the  turn  rate  is  constant).  So  the  inward  radial  body  surfaces  forces  must  increase  to 
oppose  the  increase  in  centrifugal  force.  Therefore  the  drift  angle  must  increase,  causing 
the  radial  component  of  the  propeller  thrust  to  increase  along  with  the  radial  component 
of  the  surface  pressure  forces.  So  the  location  of  the  pivot  point  will  move  closer  to  the 
nose  of  the  body. 

Wetzel  et.  al  [21]  examined  the  crossflow  separation  of  flow  over  a  submarine  in  a 
turning  manoeuvre,  and  investigated  the  effects  of  vortex  generating  fins  and  jets  on  the 

b.e  with  rounded  noses,  pointed  tails  and  a  long  parallel  mid-sections 
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location  of  this  separation  line.  Wetzel  was  limited  to  placing  the  submarine  model  in  a 
yawed  position  in  rectilinear  flow.  He  observed  changes  in  the  forces  and  moments  on  the 
model  due  the  the  delayed  separation  of  the  vortex  when  using  fins.  This  experimental 
configuration  was  able  to  produce  a  flow-field  which  simulated  a  body  in  a  gentle  turn, 
but  simply  yawing  a  body  is  unable  to  produce  the  variation  in  angle  of  incidence  along 
the  hull  section  that  is  present  when  an  object  is  performing  a  tighter  turn. 


1.6  Numerical  simulation  of  a  manoeuvring  submarine 

In  recent  years,  the  availability  of  large  computational  resources  and  developments  in 
turbulence  modelling  has  resulted  in  several  numerical  studies  of  flow  over  a  submarine 
model  (see  Kempa  et.  al  [22],  Feng  et.  al  [23],  Taylor  et.  al  [24]).  The  latest  efforts 
of  Pankajakshan  et.  al  [25],  continuing  the  work  of  Taylor  et.  al  [24]  have  resulted  in 
the  ability  to  simulate  rudder-induced  manoeuvres  of  a  submarine  model  with  a  rotating 
propulsor,  by  combining  an  unsteady  Reynolds-averaged  Navier-Stokes  equations  solver 
with  a  six  degree  of  freedom  rigid  body  model.  Their  solutions  contain  a  number  of 
interesting  features,  including  : 

1.  The  use  of  a  multiblock  dynamic  grid  system. 

2.  The  use  of  deformable  grid  to  simulate  control  surface  movement. 

3.  Incorporation  of  both  the  k  —  e  and  k  —  u  turbulence  models  on  the  hull  geometry. 

4.  The  use  of  a  body  force  propulsor  model  and  a  rotating  blade  propulsor  using  the 
k  —  ui  turbulence  model. 

1.7  The  use  of  deformed  models 

At  the  outset  of  this  project.  Professor  Peter  Joubert  suggested  a  novel  experimental 
method  which  simulates  the  flow  curvature  that  is  created  when  a  body  is  executing  a 
steady  turn  of  constant  radius.  The  experimental  model  is  itself  bent  around  the  radius 
of  turn  and  placed  in  a  straight  (rectilinear)  flow.  In  this  way,  the  variation  of  the  angle 
of  attack  (a),  along  the  length  of  the  model  is  the  same  as  for  a  turned  model.  The 
advantages  of  this  method  include  : 

1.  There  is  no  requirement  for  complicated  and  expensive  rotating  machinery. 

2.  The  object  does  not  rotate,  hence  it  cannot  pass  through  its  own  wake  (as  can  occur 
with  rotating- arm  experiments). 

3.  It  is  far  easier  to  observe  the  flow  field  and  to  record  detailed  measurements  from  a 
stationary  model  than  from  a  rotating  one. 

4.  There  is  no  time  limit  to  record  measurements,  unlike  a  rotating  arm  which  is  limited 
to  one  revolution  before  it  passes  through  its  own  wake. 
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Since  then,  a  literature  survey  conducted  by  the  authors  has  shown  that  this  method 
has  been  used  previously  to  determine  rotational  derivatives  for  aircraft,  and  to  study  flow 
curvature  effects  on  wind-turbine  blades. 

As  far  as  the  authors  can  ascertain,  the  first  recorded  use  of  this  method  was  carried 
out  at  the  EAO  TsAGI  in  Moscow  in  the  1930’s.  It  was  financed  by  the  NIO  Dirigiblestroy 
(Airship  industry)  to  investigate  the  stability  derivatives  of  airships.  The  work  was  first 
suggested  by  Prof.  Vetchinkin,  who  in  turn  refers  to  some  early  theoretical  work  conducted 
by  Glauert.  in  estimating  the  effect  of  wing  rotation  on  lift.  A  thorough  analysis  was  carried 
out  by  Gurzienko  [26]  in  1934.  It  is  also  known  that  Von  Karman2  used  this  method  to 
carry  out  experiments  on  airships  for  the  U.S.  Navy  at  the  Guggenheim  Airship  Laboratory 
in  Akron  from  1936.  This  method  was  initially  applied  to  airships  due  to  the  importance 
of  the  calculation  of  stability  coefficients  Cmq  and  Crrid  ~  the  derivative  of  pitching 
moment  (normalised)  with  respect  to  rate  of  pitch  (q)  and  rate  of  change  of  incidence 
(d).  The  small  size  of  airship  control  surfaces  means  that  the  airship  design  must  be 
inherently  stable  with  respect  to  pitching,  otherwise  the  airship  may  enter  a  pitching 
oscillation  mode  that  it  is  unable  to  correct.  Using  this  method  of  deformed  models,  these 
coefficients  are  easily  obtained.  More  work  was  done  in  Canada  during  the  1950’s  by  Foster 
and  Haynes  [27]  and  Ruckemann  [28,29]  using  these  methods  for  wind  tunnel  testing  of 
dynamic  derivatives.  As  the  use  of  more  sophisticated  wind-tunnel  apparatus  (such  as 
dynamically  oscillating  the  wind-tunnel  models)  were  developed,  these  derivatives  could 
be  obtained  without  the  use  of  curved  models.  The  method  became  less  widely  used, 
although  it  was  briefly  used  to  calculate  derivatives  of  ballistics  in  Australia  by  Forsyth 

[30]. 

A  similar  approach  to  the  problem  of  flow  curvature  effects  was  conducted  by  Migliore 
et.  al  [31],  Migliore  and  Wolfe  [32].  Migliore  investigated  flow  curvature  effects  on  Darrieus 
wind  turbines,  a  special  type  of  turbine  which  has  a  continuous  looped  airfoil  section  that 
rotates  around  its  vertical  axis.  If  the  blade  chord  of  the  airfoil  section  is  large  relative 
to  the  radius  of  rotation,  then  flow  curvature  effects  will  be  significant.  He  then  used  a 
conformal  transformation  to  transform  the  geometric  airfoil  in  curved  flow  to  an  equivalent 
virtual  airfoil  in  rectilinear  flow.  This  transformation  introduces  a  virtual  camber  and 
incidence  into  the  airfoil,  which  enabled  an  investigator  to  compare  airfoil  data  collected 
from  rectilinear  flow  experiments  and  apply  them  to  the  rotating  flow  situation.  The  use  of 
this  transformation  also  allowed  an  investigator  to  change  the  surface  profile  of  an  existing 
airfoil  section  into  a  airfoil  section  that  will  retain  the  same  lift  and  drag  properties  in  the 
rotating  flow  of  the  wind  turbine  blade.  Migliore’s  methods  were  also  examined  by  Hirsch 
and  Mandal  [33]  and  Cardona  [34]. 


2  Transforming  Equations 

Consider  a  body  performing  a  turn  around  radius  i?o,  aligned  at  an  angle  of  attack  fio 
to  the  flow  (see  figure  4).  From  triangle  BOC  we  find  : 


OB  =  Rq  cos  fio 

‘The  exact  reference  has  yet  to  be  found 


(5) 
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Figure  4:  Body  of  revolution  turning  about 
point  O.  (not  to  scale) 


Figure  5:  Body  of  revolution  deformed 
about  curved  axis.  Note  the  linear  distri¬ 
bution  of  free- stream  velocity  (deformation 
exaggerated  for  clarity) 


At  some  point  D  along  the  axis  of  the  body,  we  can  define  the  local  angle  of  attack  ft 


as  : 


tan/3  = 


BD 


xi 


(6) 


i?0  COS  fto  Ro  COS  fto 

The  model  axis  is  deformed  such  that  the  local  angles  of  attack  are  preserved  for  rectilinear 
How  (figure  5)  in  co-ordinate  system  (x,z').  In  this  co-ordinate  system  : 


dz' 

tan  ft  =  —  = 


xi 


(7) 


dx  Ro  cos  fio 

The  length  of  the  model  axis  x\  has  to  be  preserved.  Gurzhienko[26]  solved  the  equation 


for  the  curved  model  axis  by  computing  the  length  of  x\  =  BD  using  : 

~2 


x\  —  BD 


■r.  Hiy 


dx 


Combining  equation  8  into  equation  7,  then  rearranging  and  differentiating  gives  : 

d2z' 


dz'V 
dx  ) 


(8) 


(9) 


dx 2  R0 cos 

By  using  reduction  of  order,  this  can  be  solved  as  a  first-order  linear  differential  equation, 
giving  the  result  : 

z'  =  Ro  cos /80  cosh  (  —]  -  1  (10) 

L  \Rocos0oJ  J 

which  describes  the  curvature  of  the  centreline  of  the  deformed  model. 

Independent  analysis  conducted  by  Mattner  [35]  at  Melbourne  University  agrees  with 
Gurzhienko’s  methods. 
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Gurzhienko  also  derived  the  equations  which  define  the  deformed  cross-sections  passing 
through  point  D  as  : 


z  =  x  i 


cosh 


(£ 


-  XD 
Xl 


arcsinh— ) 
nx  i  / 


1 

71 


(11) 


where  n  =  Racos/3 o* 

It  can  be  shown  that  by  deforming  our  body  using  equations  10  and  11,  the  change  in 
the  local  angle  of  attack  along  our  deformed  body  is  usually  within  2-2.5%  of  our  original 
body,  despite  the  degree  of  contraction  and  expansion  along  the  inner  and  outer  surfaces 
of  the  deformed  model.  The  differences  between  the  curved  centreline  and  a  circular  arc 
segment  is  within  2%  for  values  of  fio  <  10  degrees. 


2.1  Model  construction 

In  order  to  simplify  model  construction,  it  can  be  assumed  that  the  cross-sections 
remain  planar  and  perpendicular  to  the  curved  axis.  Gurzhienko  showed  that  for  a  typical 
model,  the  radius  of  curvature  of  cross  sections  with  large  diameters  was  typically  20  times 
the  length  of  the  model.  The  advantage  of  preserving  the  circular  cross  sections  is  that  it 
makes  fabrication  of  wind  and  water  tunnel  models  far  easier.  In  addition,  creation  of  the 
deformed  surface  is  also  easier  for  use  in  numerical  simulations. 


3  Numerical  Simulations 


The  CFD  code  FLUENT  was  used  to  perform  simulations  of  flow  over  a  body  in  a 
turn  of  constant  radius.  This  was  compared  to  rectilinear  flow  over  the  corresponding 
deformed  body.  In  order  to  simulate  the  manoeuvring  body,  the  entire  flow  domain  was 
rotated  around  a  given  turn  radius.  In  this  way,  the  variation  of  magnitude  and  angle 
of  attack  of  velocity  along  the  body  could  be  replicated.  This  does  introduce  a  radial 
pressure  gradient  that  is  not  present  in  the  real  situation. 


3.1  Meshing  scheme 

The  surface  of  the  body  was  de-composed  into  logical  squares  so  that  hexahedral 
meshes  could  be  used  along  the  surface.  A  mapping  scheme  was  then  used  to  resolve  the 
boundary  layer  and  a  volume  within  half  the  diameter  of  the  model.  Beyond  this  point, 
an  unstructured  tetrahedral  mesh  was  used  to  fill  the  flow  domain. 

Since  the  boundary-layer  was  to  be  integrated  all  the  way  to  wall,  dense  mesh  spac¬ 
ing  was  required  at  the  body  surface.  Several  grids  were  generated  with  different  wall 
resolutions,  as  well  as  difference  flow  domain  sizes  to  ensure  grid-independence. 
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Figure  8:  Detail  of  surface  grid  and  boundary-layer  region 


3.2  Numerical  scheme 

Results  were  obtained  using  the  Realisable  k  —  e  turbulence  model.  The  two-layer 
zonal  method  was  used  for  wall-treatment,  which  requires  that  the  first  grid  point  in  the 
boundary  layer  to  be  in  the  viscous  sub-layer,  i.e  with  y+  <  5,  and  with  at  least  10  points 
within  Rey  <  200.  Pressure  and  momentum  were  solved  using  the  segregated  solver,  with 
second  order  accuracy. 


4  Results 

Results  have  so  far  been  obtained  for  a  body  performing  a  turn  about  a  single  radius. 
In  this  turn,  the  radius  was  three  times  the  body  length,  with  /?o  =  10.5  degrees,  and  the 
pivot  point  (point  where  /3  =  0)  was  located  at  the  bow  of  the  model.  Further  experiments 
will  examine  the  effects  of  changes  in  radius,  /?o  and  the  location  of  the  pivot  point. 

Comparisons  were  made  in  the  rectilinear  flow  cases  for  the  curved  model  between 
simulations  with  a  uniform  freestream  velocity,  and  those  with  a  shear  flow,  where  the 
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velocity  gradient  is  equal  to  : 


v  =  ujR 


(12) 


where  u  is  the  angular  velocity  of  the  model,  and  R  is  the  turn  radius  (see  figure  5).  This 
can  be  transformed  into  the  rectilinear  co-ordinate  system,  giving  : 


v  =  v0 


(13) 


For  Re  =  7.8  x  106,  the  average  static  pressure  over  the  model  is  : 


Average  static  pressure  (pa) 

Co  —  ord.  system 

Ro 

Pressure 

Rotating 

3L 

-43.6 

Rectilinear 

3L 

-42.8 

Rectilinear  (with  shear) 

3L 

-55.5 

Table  1:  Average  gauge  static  pressure  on  the  body. 


4.1  Surface  features 

Contour  plots  of  wall  shear  stress  presented  in  figure  9  along  the  rear  of  the  body  seem 
to  indicate  that  the  region  of  possible  flow  separation  (zero  wall  shear  stress)  is  larger  for 
the  curved  model. 


Figure  9:  Detail  of  wall  shear  stress  for  model  in  rectilinear  flow  (left)  and  rotating  flow 
(right) 

Plots  of  path  lines  along  the  body  show  that  flow  along  the  separation  line  of  the  curved 
body  leaves  the  surface  and  creates  a  vortex  sheet  (see  figure  10).  The  same  plot  for  the 
straight  body  in  curved  flow  (figure  11)  shows  similar  surface  patterns,  but  no  visible 
detached  vortex  sheet  is  observed.  Varying  the  velocity  magnitude  along  the  curved  body 
using  equation  13  produced  no  discernible  changes  in  the  wall  shear-stress  and  surface 
flow  patterns  when  compared  to  the  uniform  flow  case,  although  flow  separation  occured 
further  upstream  along  the  separation  line. 
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Figure  10:  Pathlines  of  particles  released  Figure  11:  Pathlines  of  particles  released 
from  curved  body  surface  in  straight  flow  from  straight  body  surface  in  curved  flow 


4.2  Flow  features 

Contour  plots  of  x-vorticity  in  figures  12,  13,  14  and  15  show  the  formation  of  the 
body  vortex  along  the  length  of  the  models  for  the  straight  body  (left),  curved  body 
(middle)  and  the  curved  body  with  velocity  gradient  (right).  Regions  of  higher  vorticity 
are  present  in  the  vortex  produced  by  the  straight  body  in  curved  flow.  Introducing  the 
velocity  gradient  across  the  curved  model  produces  better  quantitive  agreement  with  peak 
vorticity  values  when  compared  to  the  straight  model  in  curved  flow.  However,  this  causes 
a  tighter  vortex  to  form  downstream  (figure  14  and  15)  when  compared  to  the  straight 
body.  Flow  separation  and  the  formation  of  the  vortex  sheet  leaving  the  surface  can  be 
seen  in  figure  14. 


5  Discussion 


The  rotating  reference  frame  was  used  as  an  approximation  to  the  real  situation  of  a 
body  performing  a  turn  of  constant  radius. 

For  an  element  of  incompressible  fluid  in  a  rotating  reference  frame  of  constant  angular 
velocity  u;,  the  Navier-Stokes  equation  can  be  written  in  vector  form  as  : 


UU  „  _  -L  _ 

— —  +  u  •  Vu  +  2uj  x  u  = - Vp  4-  vV  u 

ot  p 


(14) 
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Figure  12:  Contours  of  vorticity  at  x  =  0.60L 


Figure  13:  Contours  of  vorticity  at  x  =  0.79L 


Figure  14 :  Contours  of  vorticity  at  x  =  0.98 L 


Figure  15:  Contours  of  vorticity  at  x  =  1.37 L 
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Figure  16:  Location  of  vorticity  measurement 


where  2a;  x  u  is  the  Coriolis  force,  and  the  centrifugal  force  u  x  (a;  x  r)  is  incorporated 
into  the  modified  pressure  p.  along  with  gravity. 

In  this  situation,  the  flow  will  have  a  pressure  gradient  of  : 


dp 

dr 


—  pruj2 


(15) 


By  using  our  co-ordinate  transformation  to  move  from  curved  flow  to  rectilinear  flow,  we 
are  changing  the  reference  frame  from  a  non-inertial  to  an  inertial  frame.  The  effects  of 
the  Coriolis  and  centrifugal  force  on  the  boundary-layer  separation  and  vortex  formation 
are  removed  when  we  make  our  co-ordinate  transformation.  The  use  of  the  curved  model 
produces  a  similar  vortex  structure  to  that  of  a  straight  model  in  curved  flow,  especially 
when  we  introduce  a  velocity  gradient  across  the  model.  However,  this  velocity  gradient 
may  not  be  able  to  be  implemented  in  an  experimental  configuration. 

The  differences  in  flow  separation  behaviour  between  the  curved  model  and  straight 
model  may  be  due  to  several  factors.  Firstly,  our  straight  model  in  rotating  flow  ex¬ 
periences  a  radial  pressure  gradient  that  would  not  be  present  for  a  manoeuvring  body 
turning  through  a  stationary  fluid.  Fluid  flowing  from  the  outboard  side  of  the  model 
to  the  inboard  side  will  experience  a  negative  pressure  gradient  which  will  help  delay 
boundary  layer  separation.  Comparing  a  model  in  rotating  flow  to  experimental  data  of 
a  manoeuvring  body  will  help  to  determine  the  effect  of  this  pressure  gradient.  Secondly, 
the  presence  of  Coriolis  and  centrifugal  forces  may  stabilise  the  flow  over  the  the  tail  of 
the  body.  These  forces  are  present  in  our  simulation  of  the  straight  model  in  rotating  flow, 
and  they  are  also  present  in  the  case  of  a  body  manoeuvring  through  a  stationary  fluid. 
Lastly,  there  may  be  Reynolds  number  effects  present  at  full  scale. 

Although  these  forces  can  be  introduced  numerically  for  a  simulation  of  straight  flow 
over  a  curved  body,  they  cannot  be  easily  incorporated  experimentally.  By  conducting 
a  numerical  simulation  of  the  curved  body  in  straight  flow,  and  including  the  Coriolis 
and  centrifugal  terms,  their  importance  in  the  validity  of  this  transformation  could  be 
ascertained. 
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6  Conclusions 


Numerical  simulations  were  conducted  comparing  a  body  placed  in  curved  streamlines 
to  a  curved  body  in  rectilinear  flow.  In  order  to  simulate  the  flow  over  of  a  body  performing 
a  turn  of  fixed  radius,  a  rotating  flow  field  was  used  to  vary  the  angle  of  attack  and  local 
velocity  along  the  length  of  the  body.  The  curved  body  was  deformed  such  that  the  local 
angles  of  attack  were  preserved.  In  this  way,  the  curved  body  was  able  to  replicate  the 
flowfield  of  the  body  in  the  rotating  fluid. 

Initial  CFD  results  show  good  agreement  between  the  two  situations.  This  suggests 
that  the  vortices  created  by  the  curved  body  will  accurately  replicate  the  vortices  created 
by  a  body  in  a  steady  turn.  Therefore,  any  experimental  results  of  the  vortices  created 
from  a  curved  model  can  be  used  with  confidence.  Some  differences  were  observed  in  the 
degree  of  flow  separation  between  the  body  in  rotating  flow,  and  the  curved  body.  These 
will  be  the  subject  of  further  investigations. 
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Appendix  A  Formulation  of  Stability  Derivatives 

Stability  characteristics  of  an  aircraft  for  example  can  be  obtained  from  an  analysis 
of  its  equation  of  motion.  These  consist  of  three  equations  describing  the  motion  of  the 
aircraft’s  centre  of  gravity  and  three  describing  the  motion  of  the  aircraft  around  its  centre 
of  gravity.  These  equations  of  motion  consist  of  mass  or  inertia  terms,  gravity,  thrust  and 
aerodynamic  terms. 

A  system  of  these  equations  usually  consists  of  three  force  and  three  moment  equations. 
Aerodynamic  forces  and  moments  are  represented  by  equations  of  motion  in  the  form  of 
series  expansions  in  variables  of  motion  and  their  time  derivatives.  The  coefficients  in 
these  expansions  are  referred  to  as  the  Stability  Derivatives. 

A.l  Example  of  hydrodynamic  stability  derivatives 

An  equation  for  the  sway  force  applied  to  a  submarine  can  be  written  as 

Y  —  m(v  —  wp  +  ur)  (Al) 

where 

v  sway  (yaw)  velocity  of  submarine 
w  heave  (pitch)  velocity  of  submarine 
u  longitudinal  velocity  of  submarine 
p  roll  rate 
r  yaw  rate 

Y  is  the  force  applied  by  the  fluid  to  the  submarine  and  includes  contributions  arising 
from  hydrodynamic  phenomena  as  well  as  from  static  stability  contributions.  Assuming 
that  it  can  be  expressed  as  a  function  of  the  instantaneous  velocities  and  accelerations  of 
the  submarine  in  six  degrees  of  freedom  and  the  deflections  of  its  control  surfaces,  then 

Y  =  Y[u,v,w,p,q,r,u,v,w,p,q,r,8b,6s,6r]  (A2) 

where 

5b  bow  plane  angle 
5S  stern  plane  angle 
5r  rudder  angle 

Using  a  Taylor  series  expansion  allows  the  fluid  forces  and  moments  to  be  given  in 
terms  of  the  sensitivity  of  each  force  and  moment  to  each  velocity,  acceleration  and  control 
deflection.  So  for  a  submarine  in  a  turn,  the  equation  for  Y  becomes 

Y  =  Yvv  +  Ypp  +  Yrr  +  Yj,v  +  Yj,p  +  Ypr  +  Ygr  8r  (A3) 
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in  addition  to  non  linear  terms  and  static  stability  terms. 
The  coefficients  are  defined  as 


(A4) 
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Appendix  B  List  of  Stability  Derivatives  in 

aerodynamic  notation 

Static  Longitudinal  Derivatives 

CmQ  pitching  moment  due  to  angle  of  attack 
Clq  lift  force  due  to  angle  of  attack 

Static  Lateral  Derivatives 

Cnp  yawing  moment  due  to  angle  of  sideslip 
Cip  rolling  moment  due  to  angle  of  sideslip 
Cyp  side  force  due  to  angle  of  sideslip 

Dynamic  Longitudinal  Derivatives 

Cmd  pitching  moment  due  to  vertical  acceleration 
Cmq  pitching  moment  due  to  pitching 
Ciql  lift  force  due  to  vertical  acceleration 
Clq  lift  force  due  to  pitching 

Dynamic  Lateral  Derivatives 

Cn/j  yawing  moment  due  to  acceleration  in  sideslip 

Cnr  yawing  moment  due  to  yawing 

Cnp  yawing  moment  due  to  rolling 

^L/3  rolling  moment  due  to  acceleration  in  sideslip 

CiT  rolling  moment  due  to  yawing 

C[p  rolling  moment  due  to  rolling 

^Y$  side  force  due  to  acceleration  in  sideslip 

Cyr  side  force  due  to  yawing 

Cyp  side  force  due  to  rolling 
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In  the  case  of  oscillatory  motion  about  a  fixed  axis,  acceleration  in  translation  (= 
rate  of  change  of  a  or  /3)  and  rotational  velocity  are  simultaneously  present.  This  gives 
oscillatory  derivatives  such  as  ( Cmq  +  Cma)  and  {Clq  +  Cl6c)  for  the  case  of  pitching 
oscillation,  and  ( Cnr  +  (Cir  +  C ja)  and  (Cyr  +  CYp)  for  yawing  oscillation. 

The  derivatives  are  usually  made  dimensionless  in  such  a  way  that: 

dCL 


etc.  where  : 


Cjn.n  — 


dCrr 


da 


.  q— >o 


Cy/3  = 


dCy 


dP 


Cnr  — 


dCn 
(j  — 

J 


CLa  = 


Cm  =  M/qSc  =  pitching  moment  coefficient 

Cl  =  L/qS  =  lift  force  coefficient 

Cn  =  N/qSb  =  yawing  moment  coefficient,  etc. 

with 

a  =  angle  of  attack 
P  =  angle  of  sideslip 
p  =  rolling  velocity 
q  =  pitching  velocity 
r  =  yawing  velocity 
V  =  air  speed 
q  —  dynamic  pressure 
S  =  reference  area 

c  =  reference  length  in  longitudinal  motion 
b  =  reference  length  in  lateral  motion 
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19.  ABSTRACT 

While  the  mechanisms  of  flow  separation  for  a  submarine  in  steady  rectilinear  flow  are  well  understood, 
the  flow  separation  behaviour  of  a  submarine  during  a  manoeuvre  is  not  as  well  understood.  The  im¬ 
portance  of  this  has  led  us  to  investigating  the  flow  separation  over  a  submarine-like  body  of  revolution 
in  a  turn  of  fixed  radius. 

The  flow  separation  over  the  manoeuvring  body  is  examined  by  using  a  deformed  body  placed  in 
rectilinear  flow  to  simulate  the  flowfield  created  by  a  body  moving  through  a  circular  path.  This  method 
was  used  by  Von  Karman  at  the  Guggenheim  Airship  Laboratory  and  by  Gurzhienko  in  Moscow  in  the 
1930s  to  investigate  the  manoeuvring  of  airships. 

The  deformed  model  is  constructed  such  that  the  local  angle  of  attack  is  the  same  for  a  regular  model 
placed  in  curved  streamlines.  Velocity  should  be  varied  linearly  across  the  test  section  to  ensure  the 
local  velocities  at  the  deformed  model  surface  are  equal  to  the  velocity  along  the  manoeuvring  model. 
Numerical  simulations  conducted  using  the  CFD  code  FLUENT  show  the  validity  of  this  transforma¬ 
tion.  Results  for  a  body  placed  in  a  rotating  reference  frame  are  compared  to  a  deformed  body  in 
rectilinear  flow. 
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